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Abstract 

Background: Fusarium wilt of banana is one of the most destructive diseases in the world. This disease has caused 
heavy losses in major banana production areas. Except for molecular breeding methods based on plant defense 
mechanisms, effective methods to control the disease are still lacking. Dynamic changes in defense mechanisms 
between susceptible, moderately resistant, and highly resistant banana and Fusarium oxysporum f. sp. cubense 
tropical race 4 (Foc4) at the protein level remain unknown. This research reports the proteomic profile of three 
banana cultivars in response to Foc4 and transcriptional levels correlated with their sequences for the design of 
disease control strategies by molecular breeding. 

Results: Thirty-eight differentially expressed proteins were identified to function in cell metabolism. Most of these 
proteins were positively regulated after Foc4 inoculation. These differentially regulated proteins were found to have 
important functions in banana defense response. Functional categories implicated that these proteins were 
associated with pathogenesis-related (PR) response; isoflavonoid, flavonoid, and anthocyanin syntheses; cell wall 
strengthening; cell polarization; reactive oxygen species production and scavenging; jasmonic acid-, abscisic 
acid-, and auxin-mediated signaling conduction; molecular chaperones; energy; and primary metabolism. By 
comparing the protein profiles of resistant and susceptible banana cultivars, many proteins showed obvious 
distinction in their defense mechanism functions. PR proteins in susceptible 'Brazil' were mainly involved in defense. 
The proteins related to PR response, cell wall strengthening and antifungal compound synthesis in moderately 
resistant 'Nongke No.T were mainly involved in defense. The proteins related to PR response, cell wall 
strengthening, and antifungal compound synthesis in highly resistant 'Yueyoukang I' were mainly involved in 
defense. 12 differentially regulated genes were selected to validate through quantitative real time PCR method. 
Quantitative RT-PCR analyses of these selected genes corroborate with their respective protein abundance after 
pathogen infection. 

Conclusions: This report is the first to use proteomic profiling to study the molecular mechanism of banana roots 
infected with Foc4. The differentially regulated proteins involved in different defense pathways are likely associated 
with different resistant levels of the three banana cultivars. 
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Background 

Musa spp. is one of the most important economic and 
agricultural crops in the world. However, the annual pro- 
duction and attribute of banana are greatly reduced by 
various infectious diseases caused by fungi, bacteria, and 
viruses. Among these diseases, the Fusarium wilt of ba- 
nana (a.k.a. panama disease) is the most important lethal 
disease; this disease is caused by the soil-borne fungus 
Fusarium oxysporum f. sp. cubense (E.F. Smith) Snyder 
and Hansen (F. oxysporum) [1,2]. The disease has caused 
severe yield losses in all banana-producing areas in Asia, 
Africa, Australia, and the tropical Americas [3]. Further- 
more, many planting areas are sharply declining because 
of F. oxysporum chlamydospores, which enable the fungi 
to persist in soil even under bad conditions and even 
in the absence of the host. Once soil is infected with 
F. oxysporum, susceptible cultivars can hardly be success- 
fully replanted in the same land for the span of 30 years 
[4]. Currently, F. oxysporum is classified into four races 
based on pathogenic characterization on different banana 
cultivars. Among the four races of F. oxysporum, race 4 
has been thought to be a major threat to the production 
of banana because it can affect almost all cultivars, except 
those infected with tropical race 1 and tropical race 2 [5]. 

F. oxysporum infection is divided into several steps: 
recognition of roots, attachment and colonization of root 
surface, penetration and colonization of root cortex, and 
hyphal proliferation of xylem vessels [6]. Particularly, the 
germination of fungal spores in soil is important in the en- 
tire process, which relies on the exudates of banana roots 
[7]. Therefore, roots are important for infection comple- 
tion and plant growth because they supply nutrients for 
fungal proliferation and assimilate water and nutrients, 
respectively. In the present study, roots were the main 
organ used for the investigations. To date, methods for 
controlling the disease include only physical and chem- 
ical measures, all of which are ineffective because the 
spread of the disease in the world has not been sup- 
pressed. Therefore, the development of resistant cultivars 
through molecular breeding based on plant defense mech- 
anisms is urgently needed, besides biological control of 
Fusarium wilt disease as one of implemented disease man- 
agement was considered [8]. 

During the course of evolution, plants have developed an 
innate immune defense system against various pathogens 
[9]. The initiation of recognition of pathogen- associated 
molecular pattern (PAMP) -triggered immunity (PTI) is 
the first branch of plant immunity, which relies on PAMP 
patterns by pattern-recognition receptors (PRRs) at the 
cell surface [10]. During the course of PTI, several intracel- 
lular responses are associated with plant defense, including 
changes in Ca 2+ flux, reactive oxygen species (ROS) and 
phytoalexin production, mitogen-activated protein kinase 
cascades, plant cell wall reinforcement at infection sites, 



and stomatal closure [11,12]. Pathogens have successfully 
evolved strategies to suppress PTI by secreting effector 
proteins directly into the plant cell to enhance virulence 
[13]. However, some of these effectors are recognized by 
plant resistance (R) proteins in the second branch of plant 
immune system to activate effector-triggered immunity 
(ETI) [14]. Generally, ETI is similar to PTI, although the 
former is believed to be stronger and faster than the latter, 
which is often accompanied by hypersensitive response/ 
localized cell death [15]. In addition, the plant immune 
system is regulated by hormone signaling molecules, such 
as salicylic acid (SA), jasmonic acid (JA), ethylene (ET), 
auxins, gibberillins, abscisic acid (ABA), cytokinins, bras- 
sinosteroids, and peptide hormones [14,16]. Among these 
molecules, SA, JA, and ET are considered important in 
plant responses to abiotic and bio tic stresses [11]. For ex- 
ample, the JA/ET pathway is involved in defense response 
to necrotrophic pathogens, herbivorous and wounding 
in Arabidopsis thaliana [17]. The ABA pathway has a 
crucial function in defense responses to drought, low 
temperature, salinity, and pathogens [18]. The auxin path- 
way has a significant function in defense responses to bac- 
teria, necrotrophic fungi, and oomycete pathogens in A. 
thaliana; in addition, the interplay of auxin and antimicro- 
bial secondary metabolites are involved in the same tryp- 
tophan pathway [19]. Antimicrobial secondary metabolites 
are mainly composed of phenolic compounds that are 
produced in the phenylpropanoid pathway to form lignins, 
lignans, neolignans, flavonoids, and anthocyanins [20]. 

Banana is an important monocot crop species, and the 
discovery of its defense mechanism is a prerequisite for 
disease control through molecular breeding. To date, re- 
search on the defense mechanism of banana against F. 
oxysporum is insufficient. Several defense-associated genes 
were identified in resistant banana roots by suppression 
subtractive hybridization [21]. Through transcriptomic pro- 
file analysis, Li et al. [22] found that the recognition of 
PAMPs and defense-related genes may contribute to Foc4 
resistance in banana. With the completion of banana gen- 
ome sequencing, proteomic profile analysis can provide a 
powerful tool to investigate the complex defense mechan- 
ism in banana [23]. In recent years, several proteomic stud- 
ies have been successfully applied to investigate the effect 
of cold tolerance and osmotic stresses on banana growth 
and development [24,25]. However, proteomic profiling of 
banana roots in response to F. oxysporum has yet to be 
performed. 

The present study investigated the resistance mechanism 
of banana roots, revealing the dynamics of defense against 
F. oxysporum infection. The defense-related genes involved 
in susceptible and resistant cultivars were identified using 
two-dimensional electrophoresis (2-DE). A total of 43 dif- 
ferentially expressed proteins in pathogen-challenged ba- 
nana roots were identified. These dynamic proteins were 
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found associated with defense. They were involved in the 
recognition of pathogen signal transduction, the produc- 
tion of ROS and antifungal compounds, the expression of 
PR genes, and the reinforcement of cell walls. This report 
is the first to use proteomic profiling to study the molecu- 
lar mechanism of banana roots infected with Foc4. This 
study revealed the dynamic changes in banana roots in- 
fected with Foc4. Our results provided insights into the 
defense mechanisms of banana and may serve as a basis 
for the control of Fusarium wilt. 

Results 

Effect of Foc4 treatment on banana roots 

In this study, we compared injury symptoms of roots in re- 
sponse to Foc4 infection in three banana cultivars (Figure 1). 
At 3 d after inoculation, no significant injury in the roots of 
highly resistant 'Yueyoukang Y was observed, and many 
progressively smaller lateral roots grew from the main 
roots. Less necrotic symptoms occurred in the roots of 
moderately resistant 'Nongke No.l', and only a few lateral 
roots grew from the main roots. However, Foc4-susceptible 
'Brazil' displayed severe necrotic symptoms in the roots. 

2-DE analysis of pathogen-induced proteins in banana roots 

To investigate the proteomic changes in banana roots 
in response to Foc4, we performed protein profile ana- 
lysis on three banana cultivars, namely, susceptible 'Brazil', 
moderately resistant 'Nongke No.T, and highly resistant 
'Yueyoukang I\ Banana roots at the fourth leaf stage were 
sampled at the time points of 3D after Foc4 or H 2 0 treat- 
ment. Subsequently, the total proteins of treated samples 
were extracted for 2-DE analysis using 18 cm gel strips. 
For each biological sample, proteomic profile analysis was 
conducted at least three times to obtain high reproducibil- 
ity. All protein spots were distributed in molecular mass 
values ranging from 10 kDa to 100 kDa, with isoelectric 
point (PI) ranging from 4 to 7. Approximately 600 protein 



spots were reproducibly matched in Coomassie Brilliant 
Blue (CBB) -stained gels. A total of 58 protein spots ex- 
hibited significant differences (at least twofold changes) in 
abundance (Figure 2). Of these protein spots, 27 were de- 
tected in susceptible 'Brazil', 16 were detected in moder- 
ately resistant 'Nongke No.l', and 15 were detected in 
highly resistant 'Yueyoukang Y (Figure 3). Particularly, 7 
separate protein spots were shared for any two cultivars, 
and 4 protein spots were shared for three cultivars (spot 
Nos. 2, 3, 7 and 13). Although most spots showed quan- 
titative changes, 5 showed qualitative changes, of which 
3 were only detected in infected susceptible 'Brazil' 
(spot Nos. 1, 11, and 12) and 2 were detected in control 
moderately resistant 'Nongke No.l' (spot Nos. 34 and 
36). A total 53 spots (91%) showed similar variation 
patterns in three banana cultivars, of which 37 were up- 
regulated and 16 were down- regulated after the infection. 
Clearly, 27 spots showed up- or down- regulation in only 
one of the three banana cultivars, of which 14 (54%), 7 
(58%), and 6 (40%) showed line-specific volume changes in 
susceptible 'Brazil', moderately resistant 'Nongke No.l', 
and highly resistant 'Yueyoukang I', respectively. These 
results implied that the proteins might be independently or 
cooperatively involved in different cellular pathways in 
response to F. oxysporum infection in banana. 

Identification of differentially expressed proteins 

To select defense proteins in banana roots that respond to 
Foc4, significant differential spots in CBB-stained gels were 
excised manually and subsequently analyzed by matrix- 
assisted laser desorption/ionization time-of-flight/time-of- 
flight (MALDI-TOF/TOF) mass spectrometer. A total of 
43 spots were successfully identified by querying the Musa 
acuminate protein database (http://banana-genome.cirad.fr/ 
content/download-dh-pahang), MASCOT NCBInr protein 
database and Plant_EST database (http://www.matrixsci 
ence.com/cgi/search_form.pl?FORMVER=2&SEARCH=SQ) 





* 1 A 







Brazil Nongke No.l Yueyoukang I Brazil NongkcNo.l Yueyoukang 1 

Figure 1 Disease symptoms from Musa roots with treatment of F. oxysporum. A: Banana roots/H 2 0 control, and B: Banana roots/infection. 
Susceptible 'Brazil', moderately resisitant 'Nongke No.l' and highly resistant 'Yueyoukang I' were treated at 3D after F. oxysporum inoculation, 
compared with water-treated control. 
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Figure 2 Representative 2-DE of proteins from Musa roots with treatment of F. oxysporum. A, Brazil/H 2 0 control; B, Brazil/infection; 
C, Nongke No.1/ H 2 0 control; D, Nongke No. 1 /infection; E, Yueyoukang I/ H 2 0 control; F, Yueyoukang l/infection. Proteins (550 \ig) were loaded 
on an 18 cm IPG strip (pH 4-7) for IEF, following electrophoresis of 12% SDS-PAGE and CBB G-250 staining. The different expression spots present 
at the position of arrows are compared with their respective controls. 



using PI and molecular weight values (Additional file 1: 
Table SI). Among these protein spots, 37 were matched to 
the M. acuminate proteins, only 1 was matched to the Jobs 
Tears (Coix lachryma-jobi L.) protein, and 5 (spot Nos. 3, 4, 
5, 6, and 7) were matched to M. acuminate expressed 
sequence tags (ESTs) without any functional information. 
Protein spot No. 35 is shown in Supporting Information Fig- 
ure 4 as an example. Most protein spots are unique. How- 
ever, more than one spot with only small shifts in molecular 
weight and/or pi were occasionally identified as the same se- 
quence in the gel, which is considered as a common feature 
of 2-DE [26]. The protein spots were identified as glutam- 
ine synthetase nodule isozyme (spot Nos. 17 and 18), puta- 
tive aconitate hydratase (spot Nos. 24 and 25), heat 
shock 22 kDa protein (spot Nos. 9 and 10), and 
Germin-like protein 12-1 (spot Nos. 11, 12, 13, and 15). 
The occurrence of protein shifts is explained as posttrans- 
lational modification or degradation. 

The predicted functions of the identified proteins were 
classified into nine categories based on the annotations 
from the Musa acuminate database, NCBInr database 
and the Gene Ontology database: defense (5%), second- 
ary metabolism (11%), polysaccharose synthesis (3%), cell 
cytoskeleton (3%), oxidative-redox stress (21%), signal con- 
duction (16%), molecular chaperones (13%), energy me- 
tabolism (3%), and primary metabolism (25%) (Figure 5, 
Additional file 1: Table SI). The protein identity of 14 up- 
regulated spots in susceptible 'Brazil' was determined: 2 
are related to disease defense reactions, 3 are related to 
oxidative-redox stress, 2 are molecular chaperone 



proteins, 1 is related to energy, 2 are involved in primary 
metabolism, and 4 are ESTs. The protein identity of 9 up- 
regulated spots in moderately resistant 'Nongke No.T was 
determined: 2 is an EST, 1 is a PR protein, 2 are related to 
oxidative- redox stress, 2 are involved in secondary me- 
tabolism, 1 is related to signal conduction, and 1 is related 
to cell polarization that is beneficial for defense. The pro- 
tein identity of 14 up-regulated spots in highly resistant 
'Yueyoukang V was determined: 4 are ESTs, 2 are related 
to defense reactions, 1 is related to cell wall strengthening, 
1 is a molecular chaperone, 1 is related to oxidative-redox 
stress, 2 are related to signal conduction, 2 are beneficial 
for primary metabolism, and 1 is involved in secondary 
metabolism. Overall, 23 up-regulated protein spots are as- 
sociated with resistance or defense in the three banana 
cultivars, except for 4 ESTs. Among them, protein spot 
Nos. 2 and 16 are for basal resistance, spot Nos. 28 and 32 
are for the production of lignin and antifungal compounds 
in secondary metabolism, spot No. 42 is for cell wall 
strengthening, spot No. 37 is for defense with cell po- 
larization, spot Nos. 13, 14, and 15 are for the ROS pro- 
duction and scavenging system, spot Nos. 33, 38, and 39 
are related to defense signal conduction, and spot Nos. 9, 
10, and 40 are molecular chaperones in various stresses. 
These results implicated that antifungal secondary metab- 
olism, lignin production, cell polarization, polysaccharose 
synthesis, and ABA defense signaling play key roles in the 
disease resistance mechanism of banana. 

The protein identity of 10 down- regulated spots in 
susceptible 'Brazil' was determined: 6 are involved in 
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Figure 3 2-DE analysis of differentially expressed protein spots 
in banana roots with F. oxysporum treatment. A: Brazil, B: Nongke 
No.l, and C: Yueyoukang I. Each sample had three biological 
replicates. Arrows indicate the abundance change of the 
identified proteins. 



to molecular chaperones, 2 are for signal conduction, and 
3 are for antioxidant defense, except for 1 EST. 

Transcript analysis of differential proteins 

In order to verify the changed expression of mRNA from 
three banana cultivars in response to Foc4 infection and 
evaluate the correlation of transcriptional and protein level, 
twelve randomly genes involved in defense, secondary me- 
tabolism, oxidative-redox stress, signal conduction, molecu- 
lar chaperones and primary metabolism were selected to 
validate through qRT-PCR method. In our study, the tran- 
script levels of the pathogenesis -related protein 1, protein 
IN2-1 homolog B, L-ascorbate peroxidase, probable gluta- 
thione S-transferase GSTF1, isoflavone reductase homolog, 
auxin-induced gene PCNT115 and 14-3-3-like gene GF14- 
6 were well consistent with their protein levels after patho- 
gen infection, while the levels of five transcripts did not 
correspond with their protein levels (Figure 6). Because 
gene expression is controlled at multiple regulatory levels at 
their transcriptional, post-transcriptional, translational, and 
post-translational levels, the validity of the correlation of 
mRNA and protein expression levels requires further study. 

Discussion 

The defense mechanism of other crops in response to 
infection has been widely studied at the molecular level. 
However, the defense mechanism of banana in response 
to infections remains obscure. To investigate dynamic 
changes in induced proteins, proteome profile analyses of 
three different banana roots were successfully performed. 
The significant differences of defense levels in banana 
roots were displayed at the time-points of 3D after Foc4 or 
H 2 0 treatment. The defense proteins acting in response to 
F oxysporum infection are involved in complex and effective 
pathways from susceptible 'Brazil' to highly resistant 
'Yueyoukang Y. A total of 43 differentially expressed proteins 
with known functions were identified by MALDI-TOF/TOF 
analysis; they were subsequently classified into nine categor- 
ies according to their functions. These induced proteins 
play important roles in defense as a response to pathogen 
infection. The elucidation of these proteins is beneficial in 
understanding the complex resistance mechanisms against 
the fungal pathogen F. oxysporum at the proteome level. 



primary metabolism, 2 are molecular chaperones, and 1 
is a signal protein, except for 1 EST. The protein identity 
of 5 down- regulated spots in moderately resistant 'Nongke 
No.l' was determined: 3 are related to oxidative-redox 
stress, 1 is a signal protein, and 1 is related to the produc- 
tion of primary metabolites. Spot 43 is the only down- 
regulated spot in highly resistant 'Yueyoukang T, and it is 
involved in primary metabolism. The protein identity of 
15 down-regulated spots in all banana cultivars was deter- 
mined: 7 is involved in primary metabolism, 2 are related 



Defense response proteins 

PR proteins are effective in responding to fungal, bacterial, 
and viral infections. Consequently, this response is related 
to the development of systemic acquired resistance (SAR). 
Proteome analyses revealed that two PR proteins involved 
in defense were induced in both resistant and susceptible 
banana after pathogen infection. The proteins were iden- 
tified as pathogenesis-related protein 1 (PR-1; No. 2) 
and RecName: Full = alpha-amylase inhibitor/endochitinase 
(No. 16). PR-1 is a dominant defense protein induced by 
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Figure 4 Identification of No. 35 protein spot from banana roots inoculated with Foc4 by MALDI-TOF/TOF MS. 



pathogens, SA, or ET, which is commonly considered as 
a marker for SAR [27]. The expression patterns of PR-1 in 
response to oomycete and fungal pathogens were posi- 
tively regulated in many plants [28]. For example, the 
antifungal protein of PR-1 in response to F. oxysporum 
was up-regulated in both GCTCV-218 and Williams [21]. 
RecName: Full = alpha-amylase inhibitor/endochitinase is a 
protein purified from the seeds of Job's Tears (C. lachryma- 
jobi L.); the antifungal protein was up-regulated to protect 
plants from insect feeding and fungal infection [29]. Two 
defense proteins were up-regulated in banana, of which the 
proteins PR-1 and RecName: Full = Alpha-amylase inhibi- 
tor/endochitinase were differentially regulated in susceptible 
'Brazil' and highly resistant 'Yueyoukang Y. Interestingly, PR- 
1 protein was only found in moderately resistant 'Nongke 
No.l\ The results indicated that more defense genes were 
expressed in response to fungal infection in highly resistant 
'Yueyoukang Y and susceptible 'Brazil' than moderately 
resistant 'Nongke No.r, where other plant defense pathways 
were regulated in moderately resistant 'Nongke No.r. 

Secondary metabolism-related proteins 

Some enzymes associated with secondary metabolism 
play key roles in the synthesis of antifungal compounds 
and cell wall components, which are involved in the phe- 
nylpropanoid pathway and other pathways. In our study, 



four proteins involved in the phenylpropanoid pathway 
were identified as caffeoyl-CoA O-methyltransferase (CC 
OMT; No. 28), isoflavone reductase homolog (so-called 
IFR "homolog"; No. 32) and leucoanthocyanidin dioxy- 
genase (LDOX; No. 34). S-adenosylmethionine synthase 
(SAM; No. 36) was identified as another protein in the sec- 
ondary metabolism pathway. The proteins in the phenyl- 
propanoid pathway were classified into three sub-pathways 
based on the following criteria: (i) isoflavonoids, (ii) flavo- 
noids and anthocyanins, and (iii) lignin biosynthesis [30]. 
The IFR "homolog" is highly similar to IFR in catalyzing 
the synthesis of isoflavonoid phytoalexin medicarpin. IFR 
in alfalfa was rapidly induced in response to fungal patho- 
gen [31]. The protein was up-regulated in highly resistant 
'Yueyoukang Y and moderately resistant 'Nongke No.l'. 
The results indicated that a large number of phytoalexins 
were synthesized to protect banana from Foc4 infection. 
CCOMT is a branch point enzyme of monolignol bio- 
synthesis, which is beneficial for the lignifications of the 
cell wall in response to pathogen infection [32]. Down- 
regulation of CCOMT in transgenic alfalfa reduces guaiacyl 
lignin biosynthesis [33]. This protein was up-regulated in 
moderately resistant 'Nongke No.r. The results implied 
that the lignifications of the cell wall have important func- 
tions in response to Foc4 infection. LDOX is a key regu- 
lated enzyme for anthocyanin synthesis, which was up- 
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Figure 5 Functional category distribution of 38 identified proteins with known functions under the treatment with Foc4. 



Li et a I. Proteome Science 2013, 11:41 
http://www.proteomesci.eom/content/1 1 /I /41 



Page 7 of 14 




Li et a I. Proteome Science 2013, 11:41 
http://www.proteomesci.eom/content/1 1 /I /41 



Page 8 of 14 



(See figure on previous page.) 

Figure 6 Transcript analysis of 12 differentially expressed genes in response to Foc4 infection by qRT-PCR. A: Pathogenesis-related 
protein 1, B: Protein IN2-1 homolog B, C: NADP-dependent malic enzyme, D: Protein disulfide-isomerase, E: Caffeoyl-CoA O-methyltransferase, 
F: L-ascorbate peroxidase, G: Probable glutathione S-transferase GSTF1, H: Superoxide dismutase [Mn] 3.1, I: Isoflavone reductase homolog, 
J: Auxin-induced protein PCNT115, K: 14-3-3-like protein GF14-6, and L: Late embryogenesis abundant protein. Each bar represents the 
mean values ± standard deviation for three biology replicates. 



regulated in response to various abiotic and biotic stresses, 
such as pathogen infection, UV light, high-intensity light, 
temperature, wounding, drought, and nutrient deficiency 
[34,35]. The protein was observed to have disappeared in 
moderately resistant 'Nongke No.l', which decreased end- 
products at the later stage of the pathway. Furthermore, the 
enzymes of the flavonoid pathway were not detected to 
change significantly in moderately resistant 'Nongke No.l\ 
Consistent with our study, the defense of moderately resist- 
ant 'Nongke No.r to Foc4 apparently does not dependent 
on the anthocyanin and flavonoid pathway. Overall, the 
different changes in the proteins of the anthocyanin and 
flavonoid pathways indicated that the plant might be 
transferring substrates between the flavonoid and antho- 
cyanin pathways to increase the production of lignin and 
isoflavonoids. SAM indirectly plays an important role in 
plant defense by catalyzing the production of S-adenosyl- 
L-methionine, which supplies as a methyl-group donor in 
the biosynthesis of numerous secondary metabolites [36]. 
This protein can increase plant resistance in response to 
various abiotic stresses [37-39]. SAM was observed to 
have disappeared in moderately resistant 'Nongke No.l\ 
These results indicated that substrates of other 
secondary metabolism pathways were transferred into 
the synthesis of lignin and isoflavonoids in response 
to Foc4 infection. 

Polysaccharose synthesis 

Protein No. 42 was identified as alpha- 1, 4-glucan-protein 
synthase (UDP-forming). UDP-forming is associated with 
the formation of cell wall polysaccharoses, including hemi- 
cellulosic and xylose, which is beneficial for the formation 
of physical barriers [40]. UDP-forming was up-regulated in 
maize in response to Trichoderma harzianum T22 [41]. Al- 
though the deposition of polysaccharoses in the cell wall 
confined the fungi to the outer root layer, the expression of 
UDP-forming was only up-regulated in highly resistant 
'Yueyoukang I'. Such results indicated that the formation of 
complex cell wall with different cell components is one of 
most effective defense strategies in response mechanisms to 
fungal degradation. This strategy was more evidently ob- 
served in highly resistant 'Yueyoukang Y than moderately 
resistant 'Nongke No.r and susceptible 'Brazil'. 

Cell cytoskeletion 

Protein No. 37 was identified as actin, which constitutes 
cell cytoskeleton that is involved in cell polarization defense 



in response to microbial attack. During the early stages of 
fungal infection, the activity of the actin pathway plays an 
important role in cell polarization defense by trafficking 
and secreting antimicrobial compounds to the infection 
site and sediment barrier material to the penetrated sites 
for cell wall enhancement, accompanied by the activity of 
defense gene PR-1 [42-44]. By contrast, in the present 
study, the protein was only up-regulated in moderately re- 
sistant 'Nongke No.l'. This result indicated that cell po- 
larization defense may be more effective for early fungal 
infection in moderately resistant 'Nongke No.l' than in 
highly resistant 'Yueyoukang Y and susceptible 'Brazil'. 

Oxidative-redox stress 

The rapid accumulation of ROS with their toxic potential, 
a phenomenon called oxidative burst, is a nearly ubiquitous 
response to abiotic and biotic stresses in plants. In the 
present study, only one protein was identified as Germin- 
like protein 12-1 (GLP, spot Nos. 11, 12, 13 and 15), which 
produces ROS. However, the high levels of ROS can lead 
to severe oxidative destruction of cell structures, such as 
nucleic acids, proteins, and lipids. To regulate oxidative- 
redox homeostasis, four antioxidant proteins were identi- 
fied as protein IN2-1 (IN2-1; protein No. 14), L-ascorbate 
peroxidase (APX; protein No.29), glutathione S-transferase 
GSTF1 (GST; protein No. 30), and superoxide dismutase 
[Mn] (SOD; protein No. 31), which were activated in 
antioxidant systems to scavenge ROS. GLP catalyzes the 
production of H 2 0 2 as oxalate oxidase or SOD in plant 
defense, which is induced by various stresses [45]. SOD 
(EC 1.15.1.1) is believed to be among the first line of 
defense proteins to scavenge ROS, which catalyzes the 
dismutation of superoxide to H 2 0 2 and oxygen [46]. The 
protein is also related to the activation of APX involved in 
the pathways of ascorbate and glutathione [47]. APX is a 
key enzyme that converts H 2 0 2 to water, in addition to 
regulating redox signal transduction [48]. Protein IN2-1 is 
identified as a homolog of GST. In addition to their func- 
tions in herbicide detoxification, stress signaling conduc- 
tion, and apoptosis regulation, GST and protein IN2-1 are 
glutathione peroxidases that function in ROS scavenging 
in various stresses [49]. The rapid accumulation of GLP 
is beneficial for resistance against rice blast and sheath 
blight [50]. The expression patterns of IN2-1/GST are 
complex under pathogen infection. For example, the in- 
fection caused by Aphanomyces euteiches (oomycota) in- 
creases the synthesis of IN2-1/GST in infected Medicago 
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truncatula [51]. By contrast, the infection caused by rust 
fungi decreases their synthesis in infected bean [52]. GLP 
were up-regulated in three banana cultivars. Meanwhile, 
IN2-1 was up-regulated in susceptible 'Brazil' and moder- 
ately resistant 'Nongke No.T. These results indicated that 
maintenance of high levels of ROS is beneficial to defense 
effectively against the pathogen infection. In our study, 
the antioxidant enzymes APX, GST, and SOD were down- 
regulated in moderately resistant 'Nongke No.T. This find- 
ing is contrary to the results of up-regulated proteins in 
response to F. oxysporum [22]. These results indicated that 
antioxidative emzyme IN2-1 were induced to defense ef- 
fectively against the pathogen infection in moderately re- 
sistant 'Nongke No.T than the antioxidant enzymes APX, 
GST, and SOD. 

Signal conduction 

Host defense mechanisms were triggered by various early 
signal pathways in response to pathogen infection. Prote- 
omic analysis showed that the signal proteins related to JA, 
ABA, and auxin in signaling pathways were differentially 
regulated in banana. Six spots were identified as signal pro- 
teins: putative horcolin (No. 1), abscisic stress-ripening 
protein 3 (ASR3; No. 8), abscisic stress ripening protein 
(ASR, No. 33), auxin-induced protein PCNT115 (AIP; No. 
35), 14-3-3-like protein GF14-E (No. 38), and 14-3-3-like 
protein GF14-6 (No. 39). Horcolin is a new jacalin-related 
lectin specific to mannose, which may play an important 
role in the perception and transduction of hormones in 
response to abiotic/biotic stress [53]. The expression of 
TaJRLLl, a jacalin-related lectin in transgenic A. thaliana, 
increases resistance to Fusarium graminearum and Botrytis 
cinema, which may be involved in the SA- and JA- 
dependent defense signal pathways [54]. ASR3 and ASR 
are stress -inducible proteins involved in ABA signal path- 
ways, which respond to various stresses, including patho- 
gen, drought, salt, cold, osmotic pressure, ABA, and injury 
[55]. ASR was up-regulated in soybean root in response to 
Fusarium solani f. sp. glycines [56]. Horcolin was induced 
in susceptible 'Brazil'. On the contrary, ARS3 were down- 
regulated in susceptible 'Brazil'. These results indicated that 
defense was regulated together by the conduction of JA 
and ABA signals. ARS was up-regulated in moderately re- 
sistant 'Nongke No.l', signifying the active defense by the 
conduction of ABA signal. AIP is related to plant defense 
and auxin-mediated plant growth. Auxin and defense- 
related antimicrobial secondary metabolites are produced 
through the same tryptophan pathway [57]. The AIP in 
grapevine was down-regulated in response to phyto- 
plasma infection [58]. In accordance with these previous 
results, the present results showed that AIP was down- 
regulated in moderately resistant 'Nongke No.l'. This 
result indicated that the substrates of the auxin pathway 
were transferred into the production of defense-related 



secondary metabolites. The 14-3-3 family comprises ubi- 
quitous proteins, which play important roles as regulators 
in plant cells; specifically, these proteins function in nitro- 
gen assimilation, sucrose synthesis, H + -pump activation, 
and ABA signal transduction [59]. Moreover, 14-3-3-like 
proteins in barley were up-regulated in response to pow- 
dery mildew fungus [60]. Accordingly, the present results 
showed that the 14-3-3-like protein GF14-E and 14-3-3- 
like protein GF14-6 were up-regulated in highly resistant 
'Yueyoukang Y. This result indicated that various plant de- 
fenses were induced upon pathogen infection. 

Molecular chaperones 

Proteins function as molecular chaperones by maintaining 
cellular homeostasis, which includes protein folding, as- 
sembly, translocation, and degradation; protein and mem- 
brane stabilization; and refolding of proteins to nonnative 
state under stressful conditions [61]. In the present study, 
four proteins were identified as molecular chaperones, 
namely, heat shock 22 kDa protein (HSP 22; protein Nos. 
9 and 10), T-complex protein 1 subunit gamma (TCP-1; 
protein No.22), protein disulfide-isomerase (PDI; protein 
No. 26), and late embryogenesis abundant protein (LEA; 
protein No. 40). HSP 22, a mitochondrial protein, is iden- 
tified as a molecular chaperone that functions in facilitat- 
ing the reactivation of denatured proteins and preventing 
the aggregation of heat-induced proteins [62]. TCP-1 aids 
in the folding and reassembly of some proteins, including 
nucleotide-binding proteins, actin, tubulin, and luciferase 
[63]. PDI acts as a molecular chaperone from the ferre- 
doxin family, whose biochemical function in plant mito- 
chondria remains unclear. However, reductive PDI has 
been associated to various roles, including removal of ab- 
errant disulfides, folding of new proteins, and reduction of 
disulfides required in activating proteins involved in anti- 
oxidant defense [64]. LEA is a hydrophilic protein that 
acts as a molecular chaperone to maintain the stabilization 
of proteins, vesicles, and endomembranes in the seques- 
tration of ions by binding water molecules [65]. In suscep- 
tible 'Brazil; HSP 22 was up-regulated, whereas TCP-1 and 
PDI were down-regulated. LEA was up-regulated in highly 
resistant 'Yueyoukang Y. These results indicated that the 
proteins had different roles as molecular chaperones in 
the defense response to pathogen infection. 

Energy metabolism 

Protein No. 21 was identified as ferredoxin-NADP reduc- 
tase. This protein catalyzes the production of NADPH to 
be supplied in C0 2 fixation; however, a substantial frac- 
tion of NADPH is involved in the electron transport of 
various cell pathways, such as N and S assimilation and 
amino acid and fatty acid metabolism [66]. Aside from 
maintaining the homeostasis of NADPH/NADP, the pro- 
tein is also considered as a key antioxidant that scavenges 
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ROS under various environmental stresses [67]. The pro- 
tein was up-regulated in susceptible 'Brazil', suggesting 
that the protein is necessary to supply energy for the 
metabolic changes in infected plant cells and to reduce 
the damage of ROS. 

Primary metabolism 

Six proteins were involved in primary metabolism: glutam- 
ine synthetase nodule isozyme (GS; Nos. 17, 18 and 19), 
putative carboxyvinyl-carboxyphosphonate phosphoryl- 
mutase (CPPM; No. 20), NADP-dependent malic enzyme 
(NADP-ME; No. 23), putative aconitate hydratase 
(Nos. 24 and 25), fructokinase-1 (FRK-1; protein No. 27), 
fructokinase-2 (FRK-2; No. 41), and malate dehydrogenase 
(MDH; No. 43). 

NADP-ME is an ubiquitous enzyme catalyzing malate 
for the production of pyruvate and NADPH, which could 
subsequently provide NADPH for the syntheses of lignin, 
flavonoid, and NADPH oxidase in plant defense [68,69]. 
NADP-ME expression is up-regulated in defense responses 
in transgenic tobacco [70]. It is rapidly induced in the 
defense response to Magnaporthe grisea [71]. The protein 
was down-regulated in susceptible 'Brazil'. Such results 
indicated that the suppression of lignin, flavonoid, and 
NADPH oxidase synthesis in susceptible banana causes 
a weak defense in response to Foc4 infection. 

GS is a key enzyme catalyzing the ammonia assimilation 
and phenyl propanoid pathway in plants, which confers 
environment stress tolerance [72]. GS was down-regulated 
in susceptible 'Brazil', indicating that the synthesis of de- 
fense compounds was increased in defense reactions. 
CPPM mainly catalyzes the formation of an unusual C-P 
bond in the processing of metabolites. To date, data linking 
this role in plants are lacking; however, the protein is well 
known for the synthesis of antibiotic bialaphos in micro- 
organisms [73]. CPPM was up-regulated in susceptible 
'Brazil'. The roles of these proteins need to be further in- 
vestigated in the future. Aconitate hydratase and MDH are 
key enzymes for citric acid metabolism, also known as 
the tricarboxylic acid cycle (TCA cycle). The protein of 
aconitate hydratase catalyzes the production of citrate 
or isocitrate [74]. MDH is also an important enzyme 
involved in multiple metabolic pathways of TCA and 
photosynthesis, catalyzing the reversible reduction of 
oxaloacetate to malate [75]. Aconitate hydratase and MDH 
were down-regulated in susceptible 'Brazil' and highly re- 
sistant 'Yueyoukang I', respectively. These results indicated 
that energy and substrates exhausted in metabolism were 
supplied by the other pathways. FRK-1 and FRK-2 are key 
enzymes of glycolysis in plants, which are both involved 
in the stress response of plant cells. Both enzymes also 
catalyze the phosphorylation of fructose to fructoses-phos- 
phate, which is a major substrate for many sugar metabolic 
pathways involved in glycolysis, starch biosynthesis, and the 



oxidative pentose pathway [76,77]. In addition, fructoki- 
nases important functions as sugar sensors or signals in 
plant defense [78]. FRK-1 was up-regulated in susceptible 
'Brazil', whereas FRK-1 and FRK-2 were up-regulated in 
highly resistant 'Yueyoukang V. These results indicated that 
the conduction of sugar signal can significantly strengthen 
plant defense. 

Conclusions 

This study is the first to perform protein profiling of 
banana roots in response to Foc4 infection. Proteome 
analyses revealed that various proteins are involved in 
complex defense pathways, including PR-related proteins, 
secondary metabolites, signal conduction proteins, cell wall 
polysaccharose synthesis proteins, cell polarization defense 
proteins, and oxidative-redox homeostasis proteins. The 
temporal expression patterns of pathogen-responsive pro- 
teins in resistant and susceptible banana cultivars may pro- 
vide useful information about the defense mechanisms 
of banana. These mechanisms were more evident in 
moderately resistant 'Nongke No.l' and highly resistant 
'Yueyoukang I' than in susceptible 'Brazil'. The proteins of 
PR and antifungal compounds were found to play main 
roles in the chemical resistance of banana cultivars against 
the fungal infection in the experiment. PR proteins were 
induced in three banana cultivars. Antifungal compound 
synthesis proteins were induced in moderately resistant 
'Nongke No.l' and highly resistant 'Yueyoukang I'. Further- 
more, other secondary metabolites were only rapidly 
induced in moderately resistant 'Nongke No.l', which was 
considered to be a more effective pathway in moderately re- 
sistant 'Nongke No.l' than in highly resistant 'Yueyoukang I'. 
Meanwhile, the structure of physical barriers was believed to 
confine the expansion of Foc4 in extracellular tissues based 
on cell wall thickness and lignifications, and polysaccharide 
precipitation. The lignin synthesis-related protein was 
down-regulated and up-regulated in susceptible 'Brazil' 
and moderately resistant 'Nongke No.l', respectively. 
Moreover, the cell wall polysaccharide synthesis-related 
protein was induced in highly resistant 'Yueyoukang I'. 
Finally, the defense mechanism is coordinated by several 
enzymes, such as signal-conductor proteins and molecular 
chaperones. Above all, PR and ROS were found to play 
main roles in defense in susceptible 'Brazil'. PR, ROS and 
antifungal compounds at the infected sites and cell wall 
lignifications play roles in defense in moderately resistant 
'Nongke No.l'. Meanwhile, PR, antifungal compounds, 
and complex cell wall proteins play roles in defense in 
highly resistant 'Yueyoukang I'. 12 randomly proteins 
involved in defense, secondary metabolism, oxidative- 
redox stress, signal conduction, molecular chaperones and 
primary metabolism were performed to validate the 
correlation of transcriptional and protein levels through 
qRT-PCR. Transcriptional levels of most selected genes 
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were consistent with their respective protein abundance 
after pathogen infection. The present study provided 
important clues for understanding the defense mechanism 
of banana against pathogen infection. Our results will be 
useful in designing strategies to cure and control diseases. 

Methods 

Plant materials and Foe 4 inoculation 

Three banana cultivars belonging to the same genomic 
group were selected in this study, including banana cv. 
'Brazil' (Musa acuminata L.; AAA Group cv. Brazil), ba- 
nana cv. 'Nongke No.T (M. acuminata L.; AAA Group 
cv. Nongke No.l), and banana cv. 'Yueyoukang Y (M 
acuminata L.; AAA Group cv. Yueyoukang I). 'Brazil' is 
susceptible to F. oxysporum, whereas 'Nongke No.l' and 
'Yueyoukang Y are moderately and highly resistant to 
this pathogen, respectively. Tropical Race 4 strain of F. 
oxysporum f. sp. cubense (Foc4) has been under long- 
term preservation in our laboratory. All cultivars were 
tissue cultured as described by Ganapathi et al. [79]. Ba- 
nana plantlets were inoculated using the method of Berg 
et al. [20] with some modifications. Propagated banana 
plantlets were transferred into 85 mm diameter tissue 
boxes containing sterile peat soil. The plantlets were in- 
cubated in a greenhouse at 21°C to 28°C with a photo- 
period of 12 h:12 h light:dark. When the plantlets 
reached the four-leaf stage, washed roots of the plantlets 
were wounded by gentle crushing prior to inoculation 
with the pathogen. Then, roots and rhizomes were 
dipped into a 500 mL sterile beaker containing 200 mL 
tropical race 4 spore suspension (10 6 conidia/mL of each 
isolate) and then incubated for 30 min. The control 
plants were treated with sterile distilled water in the 
trial. Each sample set-up was done in triplicates. Each 
treated plant was subsequently planted singly in a 200 mL 



plastic cup with steam-sterilized soil, placed on iron trays, 
and then kept in a humidity chamber in a greenhouse 
(28 ± 1°C and a natural photoperiod of 12 h:12 h light/dark). 
Root samples were collected at the point of 3D after Foc4 or 
H 2 0 treatment. The samples were flash frozen in liquid 
nitrogen and stored at -80°C. 

Protein extraction 

Proteins were extracted and separated from the samples 
according to the protocol of Isaacson et al. [80]. Briefly, 
mortar and pestle were chilled with liquid nitrogen be- 
fore grinding, and 1 g of root tissue was frozen in liquid 
nitrogen. Samples were ground to fine power with a pestle 
in liquid nitrogen and then suspended in 10 mL extraction 
buffer (0.7 M sucrose, 0.1 M KC1, 0.5 M Tris-HCl, 
50 mM EDTA and 0.2% DTT, pH 7.5). Total proteins 
were precipitated overnight at -20°C. Subsequently, 
the mixtures were centrifuged at 5,000 g for 30 min at 4°C. 
The sediment pellets were washed three or five times with 
10 mL of ice-cold methanol containing 0.2% DTT. The 
collected pellets were air-dried at 4°C and then re- 
suspended in lysis buffer [7 M urea, 2 M thiourea, 1% 
CHAPS, 2% ampholytes (pH 4 to 7), 1% DTT]. Finally, the 
mixtures were centrifuged at 16,000 g for 10 min at 4°C, 
and protein supernatants were transferred to new tubes. 
The protein concentrations were estimated using the 
Bradford assay [81]. 

2-DE and gel staining 

The lysed samples were diluted with rehydration buffer 
(7 M urea, 2 M thiourea, 4% CHAPS, 0.002% bromophenol 
blue, 0.5% IPG buffer, 75 mM DTT) to 380 uL protein 
supernatant containing 550 ug protein. Subsequently, 
18 cm IPG strips (pi 4 to 7, GE Healthcare) were 
rehydrated for at least 12 h in diluted protein. Ettan 



Table 1 Primer sequences used for qRT-PCR 



Spot no. 


Protein name 


Forward primer (5' — > 3') 


Reverse primer (5' — > 3') 


Size (bp) 


2 


Pathogenesis-related protein 1 


TOGGTACGTGAAGGAAAGA 


GGTAGGCTGATGGGTOG 


92 


14 


Protein IN2-1 homolog B 


CmOTCOTGGGCAATOAG 


TCTGGAACCTCTCAACAAATGGA 


66 


23 


NADP-dependent malic enzyme 


GCACATATOCCGCTAACG^ 


CAGCCGGGTOCCAAAC 


60 


26 


Protein disulfide-isomerase 


CCGCATCGGTGTOAGTAAA 


ACATOTCGmGCATCAACOT 


67 


28 


Caffeoyl-CoA O-methyltransferase 


GATCAACGCCAAGAACACCAT 


GGTGGCGAGGAGGGAGTAC 


61 


29 


L-ascorbate peroxidase 


GTGGACACACACTGGGAAGGT 


AGTCCAAGCCCCCTCAAAA 


58 


30 


Probable glutathione S-transferase GSTF1 


GCGAmCCGCCATCTACAT 


ACGCCGCCGAAGAAGAC 


57 


31 


Superoxide dismutase [Mn] 3.1 


TGGATAAGGGAAGTAAGAAAC 


TGCCCAGCAAAGGAAC 


93 


32 


Isoflavone reductase homolog 


CCAAGAGGCGCCAATCC 


CAAAGGCGGAGTGGCAAATA 


56 


35 


Auxin-induced protein PCNT115 


AAGGCGTGCGCATGCT 


CCACAGCGACCACTCCATCT 


58 


39 


14-3-3-like protein GF14-6 


GGATATOCCCTGGCTGAAC 


GTGCCAGCCCCAACCTAAT 


56 


40 


Late embryogenesis abundant protein 


TCGTGGATGTGCCAA^C 


GGGATCTCGCCAGTOTCTCT 


61 




Tubulin 


TGTOCATCCTGGTACTGCT 


GGCmOTGCACTGGTACAC 


112 
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IPGphor II (GE Healthcare) was used to perform IEF 
at 20°C with a current limit of 50 mA per strip. The 
program was set as follows: 12 h at 50 V, 2 h at 500 
V, 2 h at 1000 V, 7 h at 8000 V (gradient), and 50000 
Vh at 8000 V. After isoelectrical focus of proteins 
was completed, each strip was individually equilibrated for 
15 min in 8 mL equilibration solution (6 M urea, 30% 
glycerol, 2% SDS, 0.002% bromophenol blue, 50 mM 
Tris, pH 8.8) containing 1% DTT and then for 15 
min in 8 mL equilibration buffer containing 2.5% 
iodoacetamide. Second-dimension separation was carried 
out on an Ettan DALT 6 System (GE Healthcare) using 
12% SDS polyacrylamide gels. Protein spots of gels were 
visualized by CBB G-250 staining [82]. Three replicate gels 
were run for each sample. 

Image analysis 

The stained gels were scanned with ImageScanner (GE 
Healthcare), and the images were statistically analyzed using 
PDQuest software (Version 8.0, Bio-Rad, Hercules, CA, 
USA). The protein expression profiles of control tissues 
were used as a reference. Then, the matched spots were an- 
alyzed manually, and falsely matched or unmatched spots 
were corrected. Percentage spot volumes were calculated 
using Students t-test, and spots with twofold abundance 
changes were selected as significant differential proteins. 

Protein digestion 

The differential protein spots were excised from the gel. 
Tryptic digestion of the selected spots was carried out as de- 
scribed by Guo et al. [83] with some modifications. Briefly, 
each excised spot was de-stained with 100 mM NH4HCO3 
in 30% ACN. The gel pieces were minced, lyophilized, and 
then rehydrated in trypsin containing NH4HCO3 at 37°C 
overnight. After digestion, the peptides were collected, 
and the pellets were washed with 0.1% TFA in 60% (v/v) 
ACN three times to dry the remaining peptides. 

MALDI-TOF/TOF MS analysis and database searching 

The differential proteins were identified by MALDI-TOF- 
TOF as described previously [84]. The MS and MS/MS 
spectra were acquired using an ABI 4700 Proteomics 
Analyzer (Applied Biosystems, Foster City, CA, USA). 
Data search files were generated by GPS -Explorer 
Software 3.6 (Applied Biosystems), and then sequence 
similarity searches in the M. acuminate protein database 
(http://banana-genome.cirad.fr/content/download-dh-pahang), 
MASCOT NCBInr database (http://www.matrixscience. 
com/cgi/search_form.pl?FORMVER=2&SEARCH=SQ) and 
Plant_EST database were performed with MultiAgent 
Supply Chain coordination Tool (MASCOT) search 
engine 2.2 (Matrix Science, United Kingdom) with default 
parameters based on peptide mass fingerprinting PMF 
and MS/MS spectra. 



Quantitative RT-PCR analysis 

RNA samples were extracted from banana roots using E. 
Z.N.A.TM plant RNA mini Kit (Omega, USA). 0.5 ug 
total RNA as template was reverse transcription as first- 
strand cDNA using iscriptTM advanced cDNA Synthesis 
Kit (Bio-Rad) according to the manufacturers instruc- 
tions. Quantitative real time PCR was performed using 
the ssoFastTM evagreen supermix(Bio-Rad) on a Takara 
Thermal Cycler Device for Real-Time System (Takara). 
Primers listed in Table 1 are designed based on M. 
acuminate database. And a-tublin have been traditionally 
used as the endogenous reference gene [85]. 20 ul PCR 
reactions consisted of 10 ul ssoFast evagreen supermix, 2 
ul template (10 x diluted cDNA), and 20um of each 
primer (forward and reverse). Each reaction was incubated 
in PCR Strip Tubes at 95°C for 30 s, followed by 40 cycles 
of 95°C for 5 s, 58°C for 30 s. Three independent 
biological replicates were run for each sample. For 
each biological replicate, three technical replicates 
were used for qRT-PCR analyses. 

Additional file 



Additional file 1: Table SI. Identification of differentially expressed 
proteins from banana root inoculated with Foc4 by MALDI-TOF/TOF MS. 
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